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We present the general dielectric formalism for Bose- 
Einstein condensed systems in external potential at finite tem- 
peratures. On the basis of a model arising within this frame- 
work as a first approximation in an intermediate temperature 
region for large condensate we calculate the damping of low- 
energy excitations in the collisionless regime. 
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The recent experimental observation of Bose-Einstein 
condensation in alkali metal vapours [p], [@| , f§ gave a 
new impetus to the research of degenerate Bose systems. 
Much of the theoretical apparatus developed in the six- 
ties for the study of superfluid helium can be well applied 
in this case also. The present situation is even a bit sim- 
pler, as, due to the low density, interatomic collisions play 
a much less important role than in case of the helium II. 
There is, however, an additional complication due to the 
presence of the external trapping potential. 

The Bose-Einstein condensation is also a prime exam- 
ple of continuous phase transitions and spontaneous sym- 
metry breaking. Indeed, the presence of the condensate 
implies that the phase invariance of the Hamiltonian is 
absent in the ground state of the system. One of the 
most spectacular consequences of a broken continuous 
symmetry is the strong and direct coupling among dif- 
ferent correlation (response) functions. Its consequences 
have been most fully exploited by the dielectric formalism 
first introduced at zero temperature |4| , || and later gen- 
eralized to finite temperature Q , |?J (for an enlightening 
discussion of features and further references see Q). A 
most useful scheme for calculation has been developed by 
Gould and Wong || . Note that the formalism proved to 
be useful in interpreting neutron scttering experiments in 
liquid helium (see papers by H.R.Glyde and by A. Griffin 
in Jic| l and for a recent experiment [pd]| ). 

Our purpose is to generalize the dielectric formalism 
to inhomogeneous systems and apply it to trapped Bose 
condensed gases. The emphasis in this paper besides 
general considerations is put on a simple model emerg- 
ing from the general theory and the calculation of the 
damping of the lowest lying excitations in the collision- 
less regime. This is motivated by the fact that at present 
the measured damping rates of low frequency collective 
excitations JT^| are in the focus of interest, partly, 



because they are more characteristic of the Bose-Einstein 
condensation than the excitation frequencies. Recently a 
number of papers have been devoted to the calculation of 
the damping in the collisionless regime in case of homo- 
geneous systems and of trapped gases as well jl6|]- [ p2[ . 
In particular the results in ]l8fl are directly relevant to 
the experimental situation |14| . 

Let us start by summarizing the extension of the dielec- 
tric formalism to the inhomogeneous system. This for- 
malism emphasizes the importance of quantities which 
are proper and irreducible simultaneously. Having in 
mind perturbation theory and the corresponding Feyn- 
man graphs, we shall call the proper part of a quantity 
a part which cannot be cut into two by splitting a sin- 
gle interaction line. In contrast, an irreducible part will 
be a graph which cannot be cut into two by splitting a 
single line corresponding to a bare propagator (i.e., un- 
perturbed Green's function). 

As usually, we start from the Hamiltonian 



H = I d 3 rft(r) 



2m 



U(r)\ V>(r) 



+ 2 J ^ J d3 ^Hn)^Hr^)v{n,fi)^(r\)i>(r 2 ) (1) 

assuming two-body interactions. Here if)(r) is the Bose 
field operator and U (r) denotes the external poten- 
tial. Separating out the condensate part we have 
ip(r) = $o(r) + (j>(r), where $o(f) =< %j}{r) > 
stands for the condensate wave function. The an- 
gular bracket here and hereafter denotes ensemble 
averaging with the grand canonical statistical oper- 
ator e -P0-nN) i^-piH-nN) _ xhc density auto- 
correlation function is defined by x(Ti> T > ^ T ') — 



TV 



n(r, r)n(r', r') 



where n(r) = n(r) — (n(r)) = 

$ (r)ft(r) + $o(f)(/>(r) + ^(r)$(f) - ^(r)4>{r)^ is the 
density deviation operator. The field operators here 
are given in Matsubara representation, r playing the 
role of imaginary time. We shall also need the tem- 
perature Green's functions defined by Gi,i(f*, r, r',r') = 

G 2 , 2 (rW» = - (T T [<Hf,r)0t( r V; 

Gi, 2 (r,T,rV) = -(T T [0(f,r)0(rV) 

and G2,i[f, r, r', r') is obtained by replacing <j> by </>t in 
G\,%{r, r, r', r'). In the Fourier transforms G a ,p(r, r', uj) 
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and x( r l, r 2,w) the frequency oj takes on only the dis- 
crete values uo n — 2mr/(h{3) where n is an integer. The 
physically meaningful retarded correlation function and 
Green's function can be constructed by a suitable ana- 
lytic continuation [ p3[ . We shall denote both the tem- 
perature and the retarded functions by the same symbol. 

The proper part of a quantity (correlation function, 
Green's function etc.) will be denoted by a tilde over the 
corresponding symbol. The definition of the proper part 
directly implies that 

x(fi,fi,u) = x(ri,r~2,Lo) 

+j: J d 3 r 3 J d 3 r^x{r l ,rXuj)v{r A ,r 3 )x{r 3 ,f 2 ,uj) , (2) 

Eq. (||) remains valid also for the corresponding retarded 
correlation functions (i.e., after analytic continuation). 
When x(^3) i%i has a P°l e m w (and x(f 3 , r 2, w ) is non- 
singular at that u) value), 

d 3 f" 3 (sin - 7%) - i J d 3 flx{f{,fi,uj)v(fl,f" 3 ) 

xt(f 3 ) = (3) 

should hold, where ^{f^) is an eigenfunction. 

A remarkable feature of the Bose condensed system is 
the 'mixing' of correlation functions of different order, 
that is a consequence of the related symmetry breaking 
J|- This is the basis of the dielectric formalism and 
can be expressed by the generalization of one of its basic 
relationships to inhomogeneous systems as 

X(fi,r 2 ,u>) = x {r) (ri,r" 2 ,uj) 

+ J d 3 r 3 J d 3 r2A a (f{,r 3 ,uj)G a ^(fi : fl,uj) (4) 

xAp(fi,fi, uj) . 

Here x is the proper and irreducible regular part of 
the density correlation function, A Q stands for the proper 
part of the anomalous vertex A a which is due to the pres- 
ence of the condensate and represents the contribution of 
graphs with one outer particle line and one outer interac- 
tion line. G a ,p is the proper part of the Green's function. 
For repeated 'spinor' indices a summation is understood. 
Eq.([l|) allows us to derive the relation 



d 3 f 2 J d 3 i%A 1 (fi,f^,co)G ll p(f 2 ,r^,oj)Gp^ a (f3,r%oj) 
= / d 3 f2 \ d 3 r^x(n,r^,^)x^ 1 (^2,r" 3 ,Lj)A a (r 3 ,fi,uj) (5) 



Here the inverses are understood in the integral opera- 
tor sense. This equation implies the coincidence of the 
eigenvalues of the density correlation functions and the 
Green's functions. The eigenvalue equation in case of the 
Green's functions is 



- J d 3 flG a ^(f{,fl,uj) 

x ^p 0t7 (fl,fi,u) - E^rl, 7-3,0;)) ip~,{f 3 ) = <p a (f{) . (6) 

Here E^-y stands for the (irreducible) self-energy, and 
£/3, 7 for its proper part. The eigenvectors £ and ip a be- 
longing to the same eigenvalue uj are related by £(ri) = 
/ d 3 f2A a (f{,f2,uj)ip a (f2), or, the other way round, 

<p a (r~i ) = y^ 3 ^ / d 3 f" 3 J d 3 flG a .j3(fl,f2,uj) 

xAp(f^,r3 ) uj)v(ri,fl)^(fl) . (7) 

Up to now we have given expressions that are valid in 
all order of the perturbation theory. Note that in the 
Bogolyubov approximation Eq.(|(]) can be transformed 
into the familiar Bogolyubov equations for u and —v, the 
components of cp a . The same remains true in the finite 
temperature generalization of the Bogolyubov theory by 

Popov m, |. 

We turn now to present a model based upon the gen- 
eral formalism introduced above and to work out the 
damping of the low energy excitations in this framework 
in the weak-coupling regime < n > g << ksT. Here 
g = AttK 2 a/m, where a is the s-wave scattering length 
assumed to be positive and to describe the interaction be- 
tween atoms at the relevant energies. The model arising 
is the generalization to nonzero external potential that 
of treated in detail for homogeneous system in Jj| . 

The anomalous vertex A a will be approximated by 

A Q (r! , 7-2 , u) = $ (n )6(f{ -f 2 ), (8) 

where the condensate wave function < &o( r l) is assumed 
to be real. 

As for the proper Green's functions, we shall consider 
the approximation when 



(huj - H^j 



G\.2 — G2.1 
H)Gi tl {fi,f2,u)=h6(r-i- 



= 
ra). 



(9) 
(10) 



Here H = T + U(f{) - n + S$o( r D + ^gn T {f{) with 
T = — j^-A standing for the operator of the kinetic en- 
ergy and n for the chemical potential. Further, = 

/<^t(f)$(r)\ is the density of the noncondcnsate part. 

Note that 62,2(^1,^2, w) = Gi,i(r2, fl, —uj). Eq. (|io|) can 
be derived when the Hartree-Fock terms are retained in 
the self-energy graphs for the contribution of thermally 
excited particles, which corresponds to the Popov ap- 
proximation 

Inserting Eqs.(Q), (|8|) into Eq.(|3|) we get 

d 3 r^X {r) (rl,r 2 ,uj)((ri) 
-^r<Mn) J d 3 T2 Gx t i(fl,f^,u) 
+Gi,i(r3,rl, -ui)] 0o(r5)£(ra) = 0. (11) 




2 



Let us divide this equation by <j>o(fi) and apply to it the 
operator (tko — H^j {—Twj — H^j . 
We get 



2 ffj ff$o(rl)£(rl) = [TiW - H 2 



d 3 r 2 x {r) (ri,f 2 ,u,)Z(r2)) (12) 



$o(rl) ti<f>o(fi) 



The condensate wave function <&o(r) can be deter- 
mined from the Gross-Pitaevskii equation j23| 



H<P (r) = 



(13) 



Keeping only the term containing u 2 on the r.h.s. of 
Eq. (|l2"|) , neglecting x^ and using Eq.(|l3|), we arrive at 



to 



(14) 



This amounts to the hydrodynamic approximation 
known to yield accurate values for the excitation fre- 
quences when huj « [i pq |, pq |, p7] . Hereafter we 
shall be interested in this regime, and our aim is the cal- 
culation of the decay rate. This is determined in leading 
order by the term containing x in Eq. ([12"|) . Therefore, 
we consider this term a perturbation and solve Eq.jl^) 
perturbatively (still omitting the small term containing 
H 2 which does not influence the decay to this order). 

We write the frequency as u) = loq + u)\ where the cor- 
rection cji is supposed to be small compared to the lead- 
ing term luq- Similarly, the eigenfunction is decomposed 

as £(»»)= 6,(0 + . 

Collecting the corrections we get in first order 

= -u; 2 Ur) ~ -V ((Mn) f Vei(rl)) . (15) 



Multiplying Eq.(15) with £o(rI) and integrating over r~i 
the r.h.s. of the equation identically vanishes (as the 
operator appearing in the hydrodynamic equation is a 
self-adjoint one), and we are left with an equation not 
containing £i(f) any longer, i.e., with an equation for ui\ 
alone. It gives 

gu /d 8 ^/t?r3S(rl)x (r) (»^i» i 3.'^o)&(»?) na s 
^1 = ^-- r ?r=ru r^TTo • (16) 
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The damping is determined by the imaginary part of 
this expression. When evaluating Eq. ([If]) we apply the 
Thomas-Fermi approximation for the determination of 
<&o(f), as a relevant approximation in the regime where 
experiments are done. This means that in Eq.(fl3|) we 
may neglect the kinetic energy term compared to the 
others. Moreover, when solving (fl4|), the finiteness of 
the temperature will be taken into account only in the 



normalization condition of the condensate wave function. 
For a justification of this procedure see Ref. plj . Thus we 
use (Mr)f = (Mo - U(r)) /g with N /N = 1 - (T/T c f 
where N and Nq are the total and the condensate number 
of particles, respectively. 

By inserting $o (r) into Eq. ( |l4|) t he zeroth order eigen- 
function £o(f) is determined |26| , p7[ by separating 
Eq. ( pT[ ) in oblate spherical coordinates and using a poly- 
nomial Ansatz. 

The quantity x (^j r_ 2; w ) 1S approximated by the 
(dressed) loop = J2 n Ga,p(r{,r2, Un)Gp >a (f2, fi,w n - to) 
where the proper Green's function G at p(f{, f^, sat- 
isfies Eqs.(^) and (|l0|). Note that according to ( |l6| ) we 
need x^ on ly inside the condensate. In the Thomas- 
Fermi approximation applied to G one obtains for large 
condensate 



X {r) {fl,f 2 ,uj) 



(2tt) 3 



d 3 fce^ (rl - r ^x (r) (fc,w). (17) 



where j(W(t,w) stands for the bubble graph with free 
propagators investigated in [g] in detail. This means that 
the spectrum related to G is quasicontinuous in our ap- 
proximation. For the damping we need only the imagi- 
nary part of x [k,ijS), which can be calculated analyti- 
cally §. 

The expression (|l6]) for the damping was evaluated 
(using the above listed approximations) by Monte-Carlo 
integration, using 10 s points for the numerator. The 
Fourier transform has previously been evaluated numer- 
ically. The error of the evaluation of the expression ( |l6| ) 
was less than 1%. The weak coupling condition is satis- 
fied if T/T c > 0.6 taking the parameters of the JILA ex- 
periment p^ |. At temperatures close to T c the Thomas- 
Fermi approximation loses its validity. Choosing as an in- 
termediate temperature T/T c = 0.7 our result is 118 s" 1 
(115 s _1 ) for the damping rate at to = (to = 2), where 
to is the value of the L z angular momentum component, 
according to which the solutions of Eq. ( ^i|) can be classi- 
fied due to the cylindrical symmetry. This agrees within 
experimental error with the measured value fl4| and also 
with the result of pq| which was derived using approxi- 
mations different from ours. 

We have applied Eq.(||) with Eq.(||) in our calculation, 
i.e., treated the density fluctuations, which are directly 
measured. The Eqs.(|^)-(0) of the dielectric formalism 
make possible to determine the corresponding parameters 
(not presented here) of the Green's functions G a _p and 
the eigenfunctions ip a when the temperature is below the 
critical one. A basic advantage of the dielectric formalism 
is that the decoupling of the density and one-particle 
fluctuation spectra can be followed when the temperature 
increases across T c || , a fact which has been exploited in 
case of liquid helium |10|, [O. In case of trapped Bose 



gases further experiments along these lines would help 
further developing the theory. 
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